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1. Introduction 

It is generally accepted that QCD is the theory of strong interactions and thus hadrons are 
bound states of quarks and glue. Any numerical solution of QCD's bound state equations, i.e. the 
Bethe-Salpeter equation for mesons and the Faddeev equation for baryons, requires Green functions 
of quarks and gluons as input. To this end it is interesting to note that in Landau gauge the quark and 
gluon propagators are sufficiently well known for phenomenological purposes; for recent reviews 
see, e.g., [1,2]. The aim is, of course, to use these results, in particular the ones for the quark 
correlations, in the relativistic bound state equations for baryons in a similar manner as for mesons. 
However, the relativistic bound-state description for baryons is considerably less understood than 
for mesons, even on the phenomenological level. Therefore the complicated task of implementing 
our knowledge of the strong-interaction domain of QCD in a first-principle continuum description 
of baryons as Poincare-covariant three-quark bound states is only approached in small steps. 

Over the last decade one has mostly followed an approach starting from a covariant Fad- 
deev equation [3] but immediately reducing its complexity by treating the two-body T-matrix in a 
separable expansion. Due to the importance of quark-quark correlations, so-called diquarks, this 
approximative treatment has been called quark-diquark approach [4-7]. Reducing the Faddeev 
equation to an effective two-body equation allowed for the calculation of many different baryon 
observables; for a collection of recent results see, e.g., [8-12]. 

Only recently it was possible to solve the full Poincare-covariant Faddeev equation for the 
nucleon [13-15]. Its evolution with the pion mass reveals already quite interesting features, espe- 
cially when compared to corresponding lattice results. It still remains a challenging task to compute 
observables other than the mass of the nucleon. 

On the other hand, the lowest lying excited state of the nucleon, the A(1232) resonance, enjoys 
a special place in the family of baryon resonances. The A(1232) resonance has the highest produc- 
tion cross section of all nucleon resonances, and it is a spin 3 /2 particle with the same quark content 
as the nucleon. Therefore, it seems obvious to attempt also for the A(1232) the corresponding so- 
lution of the Poincare-covariant Faddeev equation. In the following we will report on progress in 
this direction by discussing the construction of the A's Faddeev amplitude and by giving a (first) 
result for its mass. 

2. Faddeev equation 

Baryons in QCD appear as poles in the three-quark scattering matrix. Analyzing the residue at 
the pole associated with a baryon of mass M allows to derive a relativistic equation for this bound 
state: 

3 

*F = tf 3 *F, K 3 =K i ; r + Y J K < $ y (2-1) 

a=l 

where now the baryon is described on-shell by the bound-state amplitude *P. This equation includes 
all possible correlations among the three quarks, which comprises a three-body irreducible contri- 
bution as well as the sum of the three possible two-body irreducible interactions with a spectator 
quark, denoted by a. 
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Figure 1: Diagrammatical representation of the Faddeev equation in rainbow-ladder truncation 



The success of the quark-diquark approach to baryon properties supports the idea that the 
quark-quark correlations dominate in binding baryons. Omitting the three-body irreducible contri- 
bution leads to the covariant Faddeev equation [3] 



a=l , 



(2.2) 



where {a,/3,y} are Dirac indices for the three valence quarks and the nature of the fourth index 
J depends on the baryon of interest. The momentum dependence of the amplitude on the three 
quark momenta p\, p 2 , pi is here reexpressed in terms of the total momentum P and two relative 
momenta p and q, where P 2 = -M 2 is fixed. They are related via: 



p = (i-OP3~tPd, q = — z — ' 

Pl=-q-2+—P, P2=q-^—P 



P-P\ +P2 + P3, 



p 3 =p + £P, 



(2.3) 



where we abbreviated Pd '■= Pi+ P2- We use the value £ = 1/3 which maximizes the upper boundary 
for the nucleon mass with respect to restrictions arising from the quark propagator's singularity 
structure; for a detailed description see, e.g., refs. [14, 16]. The internal relative momenta are given 
by: 



(1) 



p + k, 



-,(2) 



p—k, 



q-k/2, q {2) =q-kl2, q^ = q + k. 



P (3) =P, 



(3) 



(2.4) 



Furthermore, HS a ' denotes the renormalization-group invariant products of a qq kernel and two 
dressed quark propagators: 



^al'fJfJ'yy' ~ °~aa' e Kpp» 7 y» S p»p>(k 2 )S yyiki), 
A2) 

' aa'fifi'yy' 
r(3) 

^aa'fifi'yy' 



^aa'fifi'yy' ~ $ fifi'^yy" aa" S y>y(kj)S a "a'(k\), 
■^wflft'w = ^yy ,f Kaa"f}f}" S a"a'(k\)S p'fi'iki). 



(2.5) 
(2.6) 
(2.7) 



The quark propagators S involved in these expressions depend on the internal quark momenta 
ki - pi - k and kj = pj + k, where k is the gluon momentum. 
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3. Rainbow-ladder truncation 

In order to perform our calculation we have to specify the qq kernel as well as the dressed 
propagators. We require our kernel to preserve the axial-vector Ward-Takahashi identity, so that 
chiral symmetry and its dynamical breaking is correctly implemented. The simplest realization of 
such a symmetry-preserving truncation is the rainbow-ladder truncation; for its definition within 
relativistic three-body equations see refs. [14, 17]. It describes the qq kernel by a ladder dressed- 
gluon exchange 

ww^^^rt* (3-D 

which is also involved in the Dyson-Schwinger equation (for reviews see e.g., [1, 18, 19]) for the 
renormalized dressed quark propagator: 

S ~af) < yP)= Z 2{ip + m) al} + {Kaa'ppifyStfpiq). (3.2) 
Jq 

Here the bare quark mass m enters as an input and k = p- q. The color structure of the ker- 



nel generates the prefactors 2/3 and 4/3 for the integrals in Eq. ( |2.2| ) and Eq. (p.2|), respectively. 



Eq. (3T) involves an effective coupling a(k 2 ^j for which we will henceforth use the parametrization 

a(k 2 ) = mfl^je^) + a uv (k 2 ) . (3.3) 



given in [20] : 

v2 



While the first term describes the non-perturbative enhancement at small and intermediate gluon 
momenta necessary to generate dynamical chiral symmetry breaking and the dynamical genera- 
tion of a constituent-quark mass scale, the second term reproduces the one-loop QCD perturbative 
running coupling auv(k 2 ) at large gluon momentum (k 2 » JS?q CD ). The infrared behavior is param- 
eterized by the dimensionless parameter q and the infrared scale A. Hereafter we will fix the latter 
to A = 0.72 GeV which has been used in [8,20,21] to reproduce experimental results for meson 
and baryon properties. We remark that ground-state properties in this framework have turned out 
to be insensitive to the choice of q in a certain range [21]. 



4. Covariant decomposition of baryon amplitudes 

Numerical assessment of the amplitude and mass of a particular baryon requires to decompose 
the Faddeev amplitude in a suitable basis. This can be done independently of any approximation but 
based solely on first principles like Poincare covariance, spatial and time-reversal symmetries and 
charge-conjugation properties. Considering this idea, a convenient basis for the nucleon amplitude 
was obtained in Refs. [13-15]. The latter turns out to be a four-rank tensor ^a^ys {i.e., J = 5) in 
Dirac space whose dimension is 256. However, for a nucleon of positive parity and positive energy 
its amplitude can be expanded in terms of a complete orthogonal basis which only involves 64 
Dirac structures: 

64 
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Each dressing function X( in Eq. (4.1) is a dimensionless scalar factor which depends on five 
Lorentz-invariant combinations: 



p, q, Zo=PT-qr, Zi-p-P, z.2=q-P- 

In our framework a hat denotes a normalized four- vector 1 , and p^ - Tp V p y with Tp V = 5^ v - P^P y 
is a transverse projector with respect to the baryon momentum. 

It is noteworthy that the basis {O} can be divided into two sets {•SJj^Jy} whose elements 

S r ij = a®I){r i ®rj)n r (P) and P r ij = (y s ®y 5 )(r i ®T J )n r (P), (4.2) 

with 

Q r (P) = A r (P)y 5 C®A + (P), (4.3) 

involve the charge conjugation-matrix C = y 4 y 2 , the positive- and negative-energy projectors A+(P) = 
(I± p) 12, and the elements encoded in 



r, = {i, \[p!t,4],pIt,4\, 



(4.4) 



where q" t - T^Tp [v q v = T^q A T . Note that the expression above introduces a relative-momentum 
dependence in our amplitude. According to this notation, a particular basis element with index I is 
fully determined by a combination of indices i, j, r with r - ± and i,j - 1 , . . . , 4. 

In contrast to the previous case, the covariant A-amplitude is a mixed tensor in which five 
indices are involved. Four of them are spinorial whereas the remaining one is a four-vector index. 
As a consequence the amplitude can be written as 



aflyS 



(4.5) 



The spin-3/2 value of the A is then obtained by applying the Rarita-Schwinger projector 



\P) = A + (P) (t^ - i/ rT V r) , Yt = n V . 



which is restricted to positive energy and obeys the constrains P^ P^(P) = y£ (P) - 0. 

Now, in order to construct the basis for Eq. (|^), we exploit the 64 covariant structures given 
in Eq. (4.2) and proceed by multiplying S r . and f r . with all possible positive-parity combinations 



y^y 5 ®l, y 5 ®y^, 



p/y 5 1 



>p?V, 



pfy 5 ®!, 

i®pV- 



(4.6) 



>y£y 5 , y u T ®y 5 , 

This exercise generates 640 preliminary covariant structures which overcount the actual number of 
independent basis elements. Indeed, because of the following identities 

(S r ^ 



or "\ 



9 r 



(<pr \ 

ij 



{rt®y 5 )\J -(/ r r 5 ®i) * J . (W) J =(r 5 ®i) 



(s 



tpr 
'j > 



ipr 

ij > 



ij 
or 

y^ij) 



(4.7) 



We work in Euclidean space. 
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we can restrict ourselves to the upper row in Eq. (4.6) and halve the number given above to 320. A 
further reduction is achieved by contracting the former elements with I®P^ V which projects out the 
terms involving P^. Consequently, 320-64 = 256 covariant structures remain. The 64 elements 
generated from (y 5 ® y^j can be ruled out as well since they turn out to be linear combinations of 
the remaining 256 - 64 = 192 basis elements: 



TO 






( s r \ 

cy 75 ® i) '/ (i®pT) 

V ij ) 


(y 5 ®i)(^(i®^PT) 


(7 5 ®i)[^j(i®§;x v ) 



in which g = {S,f}. Only 128 of them are linearly independent. We are free to choose a linearly 
independent subset; in particular, the following choice is orthonormal: 



V5(M ~ \ . V2(M 4 + \ Wo) ■ 



(4.8) 



We note that it is rather complicated to perform a partial- wave decomposition, cf. [16], for this 
basis. A better procedure would be to apply a Gram-Schmidt algorithm and include the elements 
[3g]ij which are likely to be the dominant contributions. Such an analysis is beyond the scope of 
the present manuscript but a study in this direction is in progress. 



5. Results and outlook 

Having defined an orthonormal basis for the A-baryon amplitude, the numerical method de- 
scribed in [14, 15] can be directly applied. The solution must have the correct symmetry under the 
permutation group. Namely, both the flavour and the Poincare-Dirac parts of the amplitudes have 
to be fully symmetric under the exchange of any two quarks, as required by the Pauli principle. 
Thanks to the iterative procedure used to solve the equations, to fulfill this condition we just use 
[^g]ii> which has the correct symmetry, as the starting function. This element does not appear 



explicitly in the basis ( |4.8[ ) but can, of course, be expressed as a linear combination of the given 
basis elements. 

The resulting A mass at the physical pion mass is given in Table 1. We see that the mass 
is slightly overestimated with respect to the physical value, indicating that further structures need 
to be taken into account. Corrections to this value, though small, would arise from considering 
kernels truncated beyond rainbow-ladder and also including the 3-body irreducible contributions 
in the bound-state equation. On the other hand, the mass in this approach is only minimally below 
the quark-diquark result [10] in which only the axial- vector diquark was considered, which means 
that the diquark approximation works already very well. 
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mI dq 


^Faddeev 


1.23 


1.28 


1.26 



Table 1: A mass (in GeV) obtained from the Faddeev equation compared to the physical value and to the 
quark-diquark calculation. 

The results presented here are, however, the first piece in a more elaborated investigation. First 
of all, as said above, it would be convenient to have a partial wave decomposition of the basis to 
gain physical insight and ease comparison with quark model calculations. Also, the dependence of 
our results with the parameters in the effective coupling must be analyzed. Finally, the ingredients 
needed to calculate the A-baryon mass and structure allow us to perform the same study for the 
Q-baryon, with the only difference being the quark mass and the corresponding propagator. To 
compare with lattice data a calculation for a range of quark masses is anyhow desirable. Work in 
these directions is in progress. 
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